format 8 
factor on 

bind 
lambda   1
G   1
F   1
end


mrgp model 
0 - 2 exp(3*lambda)
0 @ 1 Erlang(50,50/4)
1 @ 3 exp(3*lambda)
2 - 4 exp(2*lambda)
2 @ 3 Erlang(50,50/4)
3 - 5 exp(2*lambda)
3 @ 0 Erlang(50,50/2)
4 - 6 exp(lambda)
4 @ 5 Erlang(50,50/4)
5 - 7 exp(lambda)
5 @ 2 Erlang(50,50/2)
6 @ 7 Erlang(50,50/4)
7 @ 4 Erlang(50,50/2)
* Reward configuration defined: 
reward
0 rew_model_0 
1 rew_model_1 
2 rew_model_2 
3 rew_model_3 
4 rew_model_4 
5 rew_model_5 
6 rew_model_6 
7 rew_model_7 
end

* Reward configuration assigned: 
bind
   rew_model_0 0
   rew_model_1 0
   rew_model_2 0
   rew_model_3 0
   rew_model_4 0
   rew_model_5 0
   rew_model_6 0
   rew_model_7 0
end


echo ********************************************************************** 
echo *********  Outputs asked for the model: model ************** 


var State_Prob prob(model, 0)
echo State probability of 0
expr State_Prob
var State_Prob_0 prob(model, 1)
echo State probability of 1
expr State_Prob_0
var State_Prob_1 prob(model, 2)
echo State probability of 2
expr State_Prob_1
var State_Prob_2 prob(model, 3)
echo State probability of 3
expr State_Prob_2
var State_Prob_3 prob(model, 4)
echo State probability of 4
expr State_Prob_3
var State_Prob_4 prob(model, 5)
echo State probability of 5
expr State_Prob_4
var State_Prob_5 prob(model, 6)
echo State probability of 6
expr State_Prob_5
var State_Prob_6 prob(model, 7)
echo State probability of 7
expr State_Prob_6


* REWARD configuration: rew
bind
   rew_model_6  0
   rew_model_5  1
   rew_model_4  1
   rew_model_3  2
   rew_model_2  2
   rew_model_1  3
   rew_model_0  3
   rew_model_7  0
end 


var Exp_SS_Reward_Rate exrss(model)
echo Expected steady-state reward rate for model
expr Exp_SS_Reward_Rate

end


